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ABSTRACT 

Context. A high percentage of the astrophysically important RR Lyrae stars show a periodic amplitude and/or phase modulation of 
their pulsation cycles. More than a century after its discovery, this "Blazhko effect" still lacks acceptable theoretical understanding. 
In one of the plausible models for explaining the phenomenon, the modulation is caused by the effects of a magnetic field. So far, the 
available observational data have not allowed us to either support nor rule out the presence of a magnetic field in RR Lyrae stars. 
Aims. We intend to determine whether RR Lyrae stars are generally characterized by the presence of a magnetic field organized on a 
large scale. 

Methods. With the help of the FORS1 instrument at the ESO VLT we performed a spectropolarimetric survey of 17 relatively bright 
southern RR Lyrae stars, both Blazhko stars and non-modulated stars, and determined their mean longitudinal magnetic field with a 
typical error bar < 30 G. 

Results. All our measurements of the mean longitudinal magnetic field resulted in null detections within 3 tr. From our data we 
can set an upper limit for the strength of the dipole component of the magnetic fields of RR Lyrae stars to ~ 130 G. Because of the 
limitations intrinsic to the diagnostic technique, we cannot exclude the presence of higher order multipolar components. 
Conclusions. The outcome of this survey clarifies that the Blazhko modulation in the pulsation of RR Lyrae stars is not correlated 
with the presence of a strong, quasi-dipolar magnetic field. 

Key words. Stars: variables: RR Lyr - Stars: magnetic fields - Polarization 



1. Introduction 

Investigations of the pulsating RR Lyrae stars contribute to al- 
most every branch of modern astronomy. These stars serve as 
standard candles to fix the cosmological distance scale and are 
often considered to be prototypes of purely radial pulsation. 
With mean periods of about half a day and brightness varia- 
tions of about one magnitude, RR Lyrae stars pulsate in the 
radial fundamental mode (type RRab), the radial first overtone 
(type RRc), and both of these radial modes simultaneously (type 
RRd). However, additional cycles occur in many RR Lyrae stars. 
A considerable fraction of RR Lyrae stars (20-30 % of Galactic 
RRab and 5% of RRc stars; |Moskalik & Poretti| |2003] |, or 
even higher numbers (about 50% according to the outcome of 
the Konkoly Blazhko survey - Jurcsik, private communication), 
shows a periodic amplitude/phase modulation on timescales of 
typically tens to hundreds of times the pulsation period. This 
phenomenon is denoted Blazhko effect (Blazhko 1907 1. More 
than a century after its discovery, the Blazhko effect remains a 
mystery. The most commonly stated hypotheses to explain the 
phenomenon invoke intrinsic effects, such as resonance, a mag- 
netic field, or variable turbulent convection. 

In the resonance models the modulation is caused by a (non- 
linear) resonance between the radial fundamental mode and a 
nonradial mode, most likely a dipole (I = 1) mode (|Cox | 1993; 
| Van Hoolst eTaT1[T998)|Dziembowski & Mizerski||2004) . Like 
the oblique pulsator model for roAp stars ( Kurtz 1982| l, the mag- 
netic models suppose that Blazhko stars have a magnetic field 
inclined to the stellar rotation axis (Cousens 1983 ; Shibahashi & 
Takat a|1995 i. The magnetic field, assumed to be a dipole field, 



deforms the radial mode leading to an additional quadrupole 



component {( = 2), for which the symmetry axis coincides with 
the magnetic axis. With the star's rotation, our view of the pulsa- 
tion components changes, causing the observed amplitude mod- 
ulation. According to this model, a magnetic dipole field with a 
strength of about 1 kG is needed for the amplitude modulation 
to be observable (Shibahashi & Takata 1995). Both the reso- 
nance and the magnetic models involve nonradial pulsation com- 
ponents, of which the presence in RR Lyrae stars still remains 
to be proven ( |Chadid et "aL]|1999t Kolenberg~| |2002|i. An alter - 
native scenario has recently been proposed by Stothers (2006). 
This model is based on purely radial pulsation and involves a 
cyclic variation of turbulent convection in the hydrogen and he- 
lium ionization zones of Blazhko stars. According to |Stofhers| 
(2006), the variation of the turbulent convection could be caused 
by a transient magnetic field in the star (turbulent dynamo mech- 
anism). The expected strength of such a field is not mentioned, 
and it would be harder to detect than a large-scale dipole field. 

Photometric observations often hint at a much more com- 
plex scenario for Blazhko stars than a modulation with one sin- 
gle period. The prototype RR Lyr shows a longer cycle of about 
4 years, beside its primary 0.567 d pulsation period, and its 40 d 
Blazhko period (which has decreased over the past decades - see 
|Kolenberg et al.|2006| l. At the end of this cycle, the Blazhko ef- 
fect suddenly weakens in strength and reappears with a change in 
the Blazhko phase ( [Detre & Szeidl|1973| >. [Szeidl| ( fT976) also re- 
ports stronger and weaker 4-year cycles in the Blazhko effect of 
RR Lyr. Other well-studied stars also reveal such longer cycles: 
about 3540 d in XZ Cyg JLaCluyze et al.||2004) , about 7200 d 
in XZ Dra Purcsik et aLl|2002) , and between 2000 and 3000 d 
in RV UMa ( Kovacs||l995| l. In addition, extensive analyses of 
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old and new photometry (same references as above) show that 
these well-studied stars exhibit coincident changes of both the 
primary pulsation period and the Blazhko period. Some Blazhko 
stars are also known to have two modulation periods, such as 
XZCyg ( |LaCluyze et1iLl|2"004"] ) and UZ UMa ( |Sodor et al~ 
2006). These photometric observations clash with the models 
requiring the Blazhko period to be exactly equal or directly pro- 
portional to the rotation period of the star, such as the latest 
version of the magnetic model (Shibahashi 2000). However, 
long-term cyclic changes could still be interpreted in terms of 
the magnetic rotator-pulsator model, by explaining the observed 
phenomena with changes of the global magnetic field structure 
and/or strength. 

From the observational point of view, the question whether 
RR Lyrae stars are magnetic is still a matter of debate. Until re- 
cently, high-precision spectropolarimetric measurements of RR 
Lyrae stars were hampered by the brightness requirements of 
the technique. RR Lyr, the prototype of the class, is by far the 
brightest with V = 7.2 - 8.2, and it is the only object of its 
kind that has been the target of spectropolarimetric observa- 
tions, the outcome of which is contradictory. |Babcock (1958) 
and |Romanov et al.| ( |1987| ) reported a variable magnetic field 
in RR Lyr with a strength up to 1.5 kG. Preston (1967]) and 
Chadid et al. (2004) reported null detections. The last group of 
authors obtained magnetic field measurements of RR Lyr with 
the MuSiCoS spectropolarimeter attached to the 2 m telescope 
at Pic du Midi (France) from 1999-2002. Having covered var- 
ious phases in the pulsation and the Blazhko cycle, they con- 
cluded from their data (<t b * 80 G) that RR Lyr is a bona fide 
non-magnetic star. 

Though no significant dipole field was detected in RR Lyr 
by |Chadid et al.| ( |2004| ), there is the need for a more exhaustive 
test of the magnetic field hypothesis to explain the Blazhko phe- 
nomenon, which can be performed through a survey of a sample 
of RR Lyrae stars. To fully strengthen or disprove the hypothesis 
of a magnetic field being the driving force behind the Blazhko 
effect, we have a strong case by also checking non-modulated 
stars. A crucial test of the magnetic models for the Blazhko ef- 
fect would be to investigate the presence of magnetic fields in 
a sample of RR Lyrae stars with different pulsational proper- 
ties. Thanks to the highly efficient instrument FORS1 attached 
to VLT unit Kueyen, numerous RR Lyrae stars have now come 
within reach for spectropolarimetry. This paper presents the first 
survey of magnetic fields in a sample of RR Lyrae stars. 

In Section 2 we discuss the targets, the observations and their 
reduction. The determination of the mean longitudinal magnetic 
field is described in Section 3. The outcome of our survey, in 
terms of the constraints on the magnetic field and the impli- 
cations for the models for the Blazhko effect, is discussed in 
Section 4, and concluding remarks are given in Section 5. 



2. Observations and data reduction 

2.1. Target selection 

We performed a high precision (cr B ^ 30 G) survey for magnetic 
fields in a sample consisting of 17 RR Lyrae stars in the south- 
ern hemisphere. The sample includes 4 non-modulated RRab 
(fundamental mode pulsators), 2 non-modulated RRc stars (first 
overtone pulsators), 10 Blazhko RRab stars, and 1 Blazhko RRc 
star. All of the targets are case studies in the understanding 
of RR Lyrae pulsation. The question we wanted to address is 
which group(s) (if any) of these stars show(s) a detectable mag- 
netic field, and in particular if a magnetic field is detected in the 



Blazhko stars (see also Kolenberg et al. 2003 ). For comparison, 
we obtained one spectrum of a magnetic A star (HD 94660). Our 
sample of stars is given in Table [T] Magnitude ranges are taken 
from the GCVS (General Catalogue of Variable Stars Kh olopov] 
|1999|), and periods from the GCVS, the GEOS RR Lyrae 



et al. 



database, LeBorgne et al. (2008]) or Kolenberg et al. (2007 
2008). 

The magnetic field strength may vary over the pulsation cy- 
cle and the Blazhko cycle, as was reported by Romanov et al. 
( |1987[ ). Based on the times of maxima given in the GEOS RR 
Lyrae database (|http : //rr-lyr . ast . obs-mip . fr/dbrr/ ), 



which contains also recent TAROT observations (ILeBorgne et al. 



2008 and references therein), and from our own data (Kolenberg 
et al.||2007| [2008), we were able to determine the pulsation 



phases of our observed stars with a precision of 0.02. For the 
determination of the Blazhko phases, in contrast, we generally 
(except for SS For) could not rely on published accurate Blazhko 
ephemerides. In order to determine the Blazhko ephemerides for 



the observed stars, we used the parameters published by Wils & 
Sodor (2005) to fit the available data of the stars, AS AS data 
as well as our own measurements (if available). Monte Carlo 
simulations were used to determine the error on the Blazhko 
period. From the available photometric data we then extrapo- 
lated the Blazhko phase at the time of our FORS 1 measurements. 
The error on the phase determination is dominated by the con- 
tribution of the error on the Blazhko period multiplied by the 
number of Blazhko cycles elapsed from the time of photometric 
measurements and of magnetic measurements. For all stars of 
our sample, the error in the determination of the Blazhko phase 
turned to be substantially higher than for pulsation phase. For 
some stars, those for which the photometric data were taken a 
few years prior to the FORS 1 measurements, the uncertainty on 
the Blazhko phase turned out to be higher than 0.5, in which case 
our phase estimate is nearly meaningless. 

Having 2.5 nights of observing time at our disposal in visitor 
mode, we optimized the observing strategy (number of different 
targets, slit width and exposure time) to get the most stringent 
constraints on the magnetic field (see Section 4.2). 



2.2. Spectropolarimetry with FORS1 

Observations were obtained in visitor mode during nights 2008- 
11-10 to 2008-11-13 with the Focal Reducer/Low Dispersion 
Spectrograph FORS 1 ( |Appenzeller et al.|1998| ) of the ESO Very 
Large Telescope (VLT). We used grism 1200 B, which covers the 
spectral range of ~ 3700 - 5100 A with a dispersion of 0.36 A 
per pixel. Slit width was set either to 0.5" or 1.0", according 
to seeing conditions, for a spectral resolution of about 2800 and 
1400, respectively. CCD readout mode was set to 200 kHz, with 
a 2x2 binning (for a pixel spatial scale of 0.25") and with a gain 
corresponding to ~ 2.2 e /ADU. 

The instrument was used in spectropolarimetric mode. 
Circular polarization measurements were performed after in- 
sterting a quarter-wave plate (at set position angles) and a 
Wollaston prism into the parallel section of the instrument's op- 
tical path. For each star, we obtained eight exposures with the 
position angle of the quarter-wave plate set at -45°, +45°, +45°, 
-45°, -45°, +45°, +45°, -45°. Such a high number of exposures 
was necessary to obtain a very high signal-to-noise ratio (SNR) 
without saturating the CCD. Setting the retarder waveplate at 
two different position angles allowed us to minimize the con- 
tributions of spurious (instrumental) polarization (e.g. Semel et 
|al.|1993[ ). The signal-to-noise ratio cumulated over both beams 
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Object Type V (mag) P(d) 



WZHya 


HIC 50073 


RRab 


10.27-11.28 





53772 


VInd 


HIC 104613 


RRab 


9.12-10.48 





47959 


ULep 


HIC 22952 


RRab 


9.84-11.11 





58148 


RUScl 


HIC 226 


RRab 


9.35-10.75 





49332 


CS Eri 


HIC 12199 


RRc 


8.75- 9.31 





31133 


MT Tel 


HIC 934/6 


RRc 


8.68- 9.28 





31690 


RVCap 


HIC 103755 


RRab-BL 


10.22-11.57 





44774 


RX Cet 


HIC 2655 


RRab-BL 


11.01-11.75 


o 


57369 


RVCet 


HIC 10491 


RRab-BL 


10.35-11.22 





62340 


RXCol 


HIC 29528 


RRab-BL 


11.4-12.6 





59404 


RY Col 


HIC 24471 


RRab-BL 


10.44-11.24 





47884 


VW Dor 


HIC 29055 


RRab-BL 


11.22-12.11 





57061 


SSFor 


HIC 9932 


RRab-BL 


9.45-10.6 





49543 


RX For 


GSC6442.00690 


RRab-BL 


11.12-12.46 





59731 


SZHya 


HIC 45292 


RRab-BL 


10.44-11.84 





53724 


XRet 


LB 3316 


RRab-BL 


11.16-12.14 





49199 


BV Aqr 


HIC 108839 


RRc-BL 


10.8-11.2 





36405 


HD 94660 


HIC 53379 


Magnetic Ap 


6.09 







and all retarder waveplate positions, estimated in the wavelength 
range 4475-4525 A, was between 1800 and 3000 per A, allow- 
ing us to reach a typical error bar of 0.05 % in the Py (the ratio 
between Stokes V and Stokes /) profile in a 1 A spectral bin. 

Data were reduced using standard IRAF routines and a ded- 
icated FORTRAN code. All the science frames were bias sub- 
tracted using a master bias obtained from a series of five frames 
taken the morning after the observations. No flat fielding pro- 
cedure was applied to the data. Spectrum extraction was per- 
formed by collapsing a 60 pixel (= 15") wide strip centred about 
the traced central peak. The extraction parameters were obtained 
from the first exposure of each series, and then adopted for all 
the remaining exposures, so as to minimize the effects of a differ- 
ent sensitivity of the CCD between the pixel regions illuminated 
by the two beams split by the Wollaston. Sky subtraction was 
performed selecting symmetric regions on the left and right side 
of each spectrum (typically between pixel 30 and 35 from the 
central peak), and fitting those with a Chebyshev polynomial. 
Owing to the limited size of the CCD region illuminated by each 
beam split by the Wollaston prism, it was not possible to calcu- 
late the sky contribution using a wider region. In fact, since our 
targets are relatively bright, sky subtraction is not a critical step, 
and could even have been skipped without significantly affecting 
the final results. 

Wavelength calibration was based on one arc frame obtained 
the morning after the observations, and typically led to a RMS 
scatter of £ 0.02 pixels. Fine-tuning of wavelength calibration 
based on night sky lines could not be performed, therefore the 
accuracy of the absolute wavelength calibration is restricted by 
instrument flexures, which are expected to be less than 1 pixel 
up to a zenith distance of 60° (see FORS1/2 User Manual). 

The final products of data reduction are the Py profiles (the 
ratio between Stokes V and Stokes /), together with their error 

l |2006) , 



bars, calculated as explained in Bagnulo et al. 



i.e., com- 



bining all eight exposures filtered using a cr-c ipping algorithm. 
In addition, we obtained the so-called "null profiles", Ny profiles 
(Donati et al. 1997), which are expected to be distributed around 
with the same FWHM that characterizes the Py profiles. A de- 



viation from of the Ny profile exceeding 3 cr would prompt for 
a re-inspection of the data and data reduction quality, and for a 
search for potential causes of spurious polarization, such as, e.g., 
those due to spectral variability during the exposures. 



It is well known that in RR Lyrae envelopes, especially those 
of RRab type stars, strong acoustic waves occur, and shock 



waves are produced at certain pulsational phases (Preston, Smak 
|& Paczynski [1963] |Chadid et al. 12008} and references therein 
for recent observations). The gas dynamics associated with these 
waves in the line forming region, as well as the existence of the 
preheating zones ahead of the shock fronts and the cooling zones 
behind them, can strongly affect the shape of the spectral lines 
in RR Lyrae stars, resulting in line asymmetry, additional broad- 
ening, line profile doubling, and emission components. These 
distortion effects are particularly present during the rise to max- 
imum light, around the phase of minimum radius. 



As pointed out by |Wade et al~| ( |2002| i and |Chadid et aL 



(2004 1, the rapid changes in the line profiles of both the Balmer 
and the metallic lines may lead to spurious polarization sig- 
nals in high-resolution spectropolarimetric data. However, it is 
reasonable to assume that this problem is less critical with the 
technique based on low-resolution spectra adopted in this work. 
Indeed, at a spectral resolution of 1400 or 2800, as was obtained 
with FORS 1 for this study, most metal lines are not resolved, let 
alone their variations. Most of our observations were obtained 
outside of the phase range in which the stars are on the steep 
rising branch of their light curve (typically the pulsation phase 
interval (f> = 0.9 - 1.0). For the few stars in our sample of which 
we obtained spectra close to the shock phase, we could not see 
any effects on the spectra, in particular no spikes show up in 
the null profiles, as we would expect if the resolved line profiles 
were significantly changing during the observing time. 
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HD 94660 MT Tel SS For 




WAVELENGTH (A) WAVELENGTH (A) WAVELENGTH (A) 



<B z ) = -1964±33 G (fl a } = -44±39 G <B 2 ) = -5±20 G 




-4.67 10" 13 X s (1/1) (df/dX) (G _1 ) 



Fig. 1. Spectropolarimetric observations of the magnetic star HD 94669 (left panels) and two program RR Lyrae stars (MT Tel, 
middle panels, and SS For, right panel). The top panels show the spectral region around Hy. The top, solid lines represent Stokes /, 
expressed in arbitrary units (and not normalized to the continuum); middle lines show the Py profiles, and the bottom lines show 
the null profile Ny, which represents an overal estimate of the noise of our circular polarization measurements. The bottom panels 
show the corresponding best fits to the full spectrum obtaining by minimizing the^f 2 of Eq. J3J. 



3. Determinations of the mean longitudinal 
magnetic field 

Our determinations of the mean longitudinal field (B z ) i.e., the 
component of the magnetic field along the line of sight, averaged 
over the stellar disk, are based on the use of the formula 

- = -g eff C, A 2 - — (B z ) , (1) 
/ 5 z I dA z 

where g e g is the effective Lande factor and 

e 



C, = 



4mn e c 2 



(a 4.67 x lO^A" 1 G' 1 ), 



(2) 



where e is the electron charge, m e the electron mass, and c the 
speed of light. 

A least-squares technique was used to derive the longitudinal 
field via Eq. ([T]), by minimizing the expression 

(yi-<B z )xt-b? 



X' 



-Z 



(3) 



where, for each spectral point i, y, = (V//),, x, = -g e ffC z A 2 {\ 1 1 x 
dI/dA)i, and b is a constant. 



Note that Eq. ([TJ is strictly correct under the so-called weak 
field approximation, which, e.g., in the atmosphere of main se- 
quence A-type stars, is valid for field strengths S5 1 kG for metal 
lines, and up to ~ 10 kG for H lines. 

Using low resolution spectra such as those obtained in this 
work, it would make sense to hypothesize that only broad H 
Balmer lines could be used to determine the magnetic field, 
since most metal lines are not resolved in the spectra. However, 
|Bagnulo et~aL] < |2002[ > and |Bagnulo et al.| ( |2006) have shown that 
Eq. ([ljl may also be applied to spectral regions with non resolved 
metal lines, leading to results consistent to the (B z ) determina- 
tions obtained on H Balmer lines in the weak field regime. We 
checked that with our spectra the (B z ) determination obtained 
from metal lines only are consistent, within the error bars, with 
the (B z ) values obtained from Balmer lines only. For Balmer 
lines we adopted an effective Lande factor of 1 (Casini & Landi 



Degl Tnnocenti| 1994 1; for metal lines, following Chadid et al. 
(2004), we used an effective Lande factor of 1 .38. Finally, we ap- 
plied the least-square technique to the full spectrum, i.e., using 
both metal and Balmer lines. This method allowed us to bring 
the (B z ) error bars down to a typical level of 20-30 G. 
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Similarly to the application by Ba gnulo et al. ( 2006| l, the (B z ) 
values were determined using the Py profile obtained from the 
combinations of all exposures. For the purpose of checking our 
results, we also calculated the average of the (B z ) values ob- 
tained from individual pairs of exposures for each star, always 
finding results consistent with the {B : ) value obtained from the 
combined Py and / profiles. 

Since FORS 1 is not equipped with circular polarization fil- 
ters that allow one to routinely check the correct alignment of 
the polarimetric optics, and no standard stars for circular po- 
larization are available, we decided to perform a health check 
by observing the well known magnetic Ap star HD 94660 
(=KQ Vel=HIC 53379) during our observing run. Two con- 
secutive series of eight observations were obtained, one with a 
0.5" slit width, and one with a 1.0" slit width, using the same 
strategy adopted for the program stars. HD 94660 presents a 
roughly constant longitudinal magnetic field, and the measured 
value of ~ -1960 G (from Balmer lines only) is in reasonably 
good agreement wit h previous observations (see, e.g., Fig. 3 of 
Bagn uhTet al.|2 006, the new observation presented in this paper 



would be located at phase 0.96.). 

Figure [T] shows our observations of stellar spectra in circu- 
lar polarization, and their interpretation in terms of magnetic 
field for HD 94660, a known magnetic star, and two represen- 
tative RR Lyrae stars: MT Tel, an RRc star, and SS For, an RRab 
Blazhko star. Note that the hotter RRc star clearly has less metal 
lines in its spectrum than the RRab star. The bottom panels show 
the best fits to the full spectrum obtained through minimizing 
the x 1 °f Eq. Whereas a clear slope is observed in the case 
of HD 94660, corresponding to the mean longitudinal magnetic 
field (B z ) , for the RR Lyrae stars the best fit results are consistent 
with a curve with no significant slope. 

The results of our survey are reported in Table [2] Inspection 
of this table show that in all cases our (B z ) measurements are, 
within 3cr, fully consistent with a null detection of the mean lon- 
gitudinal magnetic field. In other words, for the 17 stars in our 
sample and the 20 measurements taken, we obtained no field 
that differed from zero exceeding the 3 <x significance level. In 
fact, 14 of our 20 measurements were null detections within a 
lcr significance level. A histogram reflecting our measurements 
is given in Fig. [2] 

4. Discussion 

4. ) . The magnetic model challenged 

The presence of a strong dipole magnetic field is a key ingredi- 



ent for the model proposed by |Shibahashi & Takata ( 1995 i and 
Shibahashi (2000). Assuming that RR Lyrae stars have fairly 
strong dipole fields, |ShibahashT ( 2000[ > showed that the radial 
mode in RR Lyrae stars would be deformed by the Lorenz force 
to have additional quadrupole (( = 2) pulsation components. 
The changing aspect angle as a consequence of the rotation of a 
star would then give rise to an apparent amplitude variation. The 
quadrupole components induced by the magnetic field would 
give rise to a quintuplet structure in the frequency spectrum. For 
more than a decade after the magnetic model was first presented, 
no evidence for quintuplet components was found from any data 
set. Shibahashi & Takata ( 1995 1 showed that a quintuplet struc- 



ture may manifest itself as only a triplet depending on the ge- 
ometrical configuration (angles of pulsation axis, magnetic axis 
and aspect angle). Also, the quintuplet components might have 
such a low amplitude that they are hidden in the noise of the 
frequency spectrum of time series data. The latter explanation 



was supported by recent findings from high-quality data sets by 
Hurta et al. (2008) in RV UMa , Jurcsik et al. (2008) in MW 
Lyr, and [Kolenberg et al. ( 2008[ ) in SS For, one of the stars in 
our sample. However, the detection of quintuplet structures does 
not imply that the magnetic model is the one to be preferred 
over the resonance model. Jurcsik et al. ( 2008a| l and Kolenberg 
|et al.| ( |2008| l also found evidence of even higher order multiplet 
components (septuplets, nonuplets) in the frequency spectra of 
Blazhko stars. These frequency structures have never explicitly 
been considered in any of the models proposed to explain the 
Blazhko effect. 

Our sample of stars consisted of RR Lyrae stars with differ- 
ent pulsational properties: RRab, RRc stars, and modulated star 
of both types. In none of these stars a field exceeding a signifi- 
cance level of 3<x was detected. 

For the magnetic field to have the effects described by 
|Shibahashi & Takata| ( [T995] l and |Shibahashi ] (|2000), the pres- 
ence of a dipole field with a strength of the order of 1 kG was as- 
sumed. In the following section we discuss what constraints our 
measurements set to the magnetic field strength and morphology 
of our sample of stars. 

4.2. Constraints on the magnetic field 




-5 -4 -3 -2 



-10 12 3 



Fig. 2. Distribution of the (B z ) values normalized to their error 
bars for our sample of 20 field measurements in 17 RR Lyrae 
stars. The dotted vertical lines correspond to the 3-<x detection 
level. 



We consider a magnetic configuration produced either by a 
dipole, or by a (non-linear) quadrupole, or by the superposition 
of both. The model is assumed to be centered, i.e. the elementary 
multipoles that produce the magnetic configuration at the stellar 
surface are located at the star's center. The dipole field Sd{r) can 
be characterized by the dipole field strength B^ at the pole, plus 
the direction of a unit vector u; the quadrupole field S q (r) can be 
characterized by the quadupole field strength Z? q , plus the direc- 
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Table 2. Determination of the mean longitudinal magnetic field of 17 RR Lyrae stars. 



Object 


Tvnp 


Date 


Time (UT) (I) 


Exp {1) 


Slit 






(B ) 

\ lJ z/ 


(SNR) CT 






fvvvv mm dd) 
\j j j j 111111 uu ; 


(hh;mm) 


(sec) 


width 


(±0.02) 




(G) 




WZHya 
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(5) The signal-to-noise ratio (SNR) is calculated per A, in the wavelength interval 4475-4525 A. 



tions of two unit vectors U\ and u^. The geometric configuration 
of a dipole plus quadrupole magnetic field is illustrated in Fig. 1 
of |Landolfi et al. (1998). Denoting by R* the star's radius, the 
magnetic field vector at a given point r of the stellar surface is 
given, respectively, by 

fl d (r) = -f [«-3^r] 

»qW = -I hr "i + 



Ui 



(4) 



( OLE 



A" 



l)r]. 



Note that while the quantity B&, the dipole field modulus at the 
magnetic pole, has an obvious physical meaning, the quantity 
B q represents the quadrupole field modulus at the magnetic pole 
only when the unit vectors U\ and u*l coincide (linear quadrupo- 
le), otherwise the concept itself of magnetic pole becomes mean- 
ingless, and the estimate of the actual field strength requires an 
evaluation of Eq. Q for the specific quadrupole configuration. 
For a more detailed description of the multipolar expansion of 
the magnetic field see Bagnulo et al. ( 1996 1. 
Assuming a limb-darkening law of the form 



1 -u + ufx (0 < u < 1) 



where ft is the cosine of the angle between the line of sight 
through the star disk center and a given point of the stellar disk, 
and using Eqs. (24), (28) and (68) of Bagn ulo et al.| ( [T996| l, we 
find that the mean longitudinal field is given by 



(B z ) = 



15 + » j 
60 1 



+ f2^q( cos cos ^2 + 3 sin £\ sin £% cos(Ai - A2)) 



(6) 



where t is the angle between the line of sight and the dipole 
axis, t\ and {2 are the angles between line of sight and the unit 
vectors u\ and «2 that define the quadrupolar component of the 
magnetic field, and Ai and A2 are the azimuth angles of the same 
unit vectors, reckoned counterclockwise in the plane of the sky 



from the north celestial meridian, as detailed in Bagnulo et al. 
fll996l 

Even a quick glance at the results in Table|2]tells us that, if a 
dipole field is present, its strength must be small, much smaller 
than the 1 kG field that the magnetic model ( Shibahashi 2000 ) 
requires to be responsible for the Blazhko effect. Based on our 
measurements, it is possible to obtain more precise constraints 
on the magnetic field of our program stars. 

Let us consider the case of a pure dipolar field. We assume 
that all stars of our sample have the same dipolar strength, and 
that the magnetic axis is randomly oriented with respect to the 
line of sight. The probability P{() for the tilt angle of the field I 
to be in the range \t, I + d{] is proportional to sin td£. A simple 
application of Bayesian statistics allows us to evaluate the prob- 
ability !P(fld ^ B^) that all the observed stars of our sample 
are characterized by a dipolar field strength higher than a certain 



threshold value B 



(0). 



HB d > Bf>) (x 



>(0> 7=1 1 



^(B ; )j-k(u)B,, cost) 



(7) 



where the product runs on the N - 11 observed stars (for the (B z ) 
value of the two stars that were observed more than once, we 
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consider the value measured first), crj is the error bar associated 
to the measurement of the / h star, and we have set 



k(u) 



15 + u 
20 (3 -it) 



(8) 



Using a hmb-darkening coefficient u = 0.62, derived from|Diaz- 



Cordoves et al. | ( [1995| ) (see also [Kolenberg 



2002), we obtain 

that if the RR Lyrae stars of our sample are all characterized by 
the same dipolar strength, with magnetic axis randomly oriented 
about the line of sight, then there is a 95 % probability that their 
dipolar strength is % 130G. 

Equation (|6]i shows that, for the limb-darkening coefficient 
of u — 0.62, the contribution of the quadrupolar component to 
the mean longitudinal magnetic field is always ^ 20 % the con- 
tribution due to the dipolar component (assuming = B q ). For 
a linear quadrupole, i.e., with u\ || «2, we always have l\ = 1% 
and Ai = A2, and our observations rule out a quadrupolar field 
strength ^ 250 G at the 95 % confidence level. If we consider 
the special case where I2 - t\ + 90° and Ai = A2 (for which 
u\ ± 112), we obtain the value B q - 1.9 kG as an upper limit to 
the quadrupole strength. 

Therefore, determinations of the mean longitudinal field as 
obtained in this work set fairly loose constraints on the magnetic 
field strength for morphologies of higher order than the dipole 
field. 



4.3. Variability over pulsation cycle and Blazhko cycle 
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Fig. 3. Absolute value of the mean longitudinal field versus pul- 
sation phase. Different symbols were used for different types of 
RR Lyrae stars. Data points for stars observed more than once 
are denoted with filled symbols. 



Figure [3] shows the plot of the absolute value of the mean 
longitudinal field as a function of the pulsation phase. Although 
this plot refers to different stars with different geometric con- 
figurations, it could potentially reveal the presence of a cyclic 



variability of the magnetic field over the pulsation cycle as was 
also observed on several nights by Romanov et al.| (|T987 ) (their 
Figs. 1 and 2). Indeed, Figure [3] shows a hint for higher (B z ) val- 
ues (though all consistent with zero) close to phase zero. This 
could be explained by the presence of a small magnetic field 
and assuming magnetic flux conservation. (The stellar radius and 
surface is minimum shortly before phase zero, and attains its 
highest value around phase <p = 0.5.) Further conclusions may 
be drawn only after a detailed monitoring of individual stars. 

For two of the stars in our sample, CS Eri (RRc) and SS 
For (RRab-BL), we obtained more than one measurement at dif- 
ferent phases in their pulsation cycle. No clear variation with 
the pulsation phase exceeding the error bars could be derived 
for these stars. We observed SS For at 3 different phases: one at 
maximum light, a second one approaching minimum light, and a 
third one on the rising branch of the light variation. From the last 
spectrum a higher absolute (B z ) was derived, but also a higher er- 
ror bar, which can be attributed to the lower signal-to-noise ratio 
of the spectrum. SS For shows quintuplets in its pulsation spec- 
trum ( Kolenberg et al. 2008 ) which, according to the magnetic 
model (Shibahashi 20001 are associated with the effects of a 
magnetic field (see also Section 4.1), and hence was one of our 
prime targets for investigating the presence of a magnetic field. 

|Romanov et al.| ( |1987| l also reported a variation of their de- 
rived mean longitudinal magnetic field values over the Blazhko 
cycle of RR Lyr. Chadid~et al.| (|2004]), in contrast, with 27 mea- 



surements of the same star, detected no significant changes at 
different Blazhko phases. Within the allocated telescope time we 
were not able to follow one Blazhko star at different Blazhko 
phases. Only for some stars in our sample the Blazhko phase 
could be determined with a reasonable accuracy. However, we 
can safely say that the 11 Blazhko stars in our sample were 
recorded at random Blazhko phases, and for all of them we ob- 
tained a mean longitudinal magnetic field consistent with a null 
detection within 3cr. As discussed in Section 4.2, this implies 
low values for the magnetic dipole field. 

Our observations significantly enlarge the sample of spec- 
tropolarimetric measurements of RR Lyrae stars and follow the 
line of results obtained by Preston (1967) and Chadid et al. 



(2004). The latter authors declared RR Lyr to be a bona fide non 
magnetic star on the basis of high-resolution spectropolarimetric 
data obtained over a 4-year time span. 

Now why did earlier measurements of the magnetic field 
in specifically RR Lyr, the brightest star of the class, lead to 
reported values of the longitudinal magnetic field strength as 
high as about 1.5 kG ( |Babcock|[T958[ |Romanov et al.||1987|), 
while more recent observations contradict this? IChadid et "aLl 
(2004) addressed this question in the necessary detail (in their 
Section 3.1), pointing at the limitations of the use of photo- 
graphic plates resulting in underestimated error bars, and the 
distortions of the spectral lines caused by shock waves in cer- 
tain pulsation phases (see also Wade et al. |2002|). The measure 



ment with by far the highest (B z ) detection by |Romanov et al. 



( |1987[ l was recorded at the phase of the main shock in the star 
(0=0.95). This may have led to spurious polarization signals. We 
also note two limiting factors of those earlier measurements: i) 
a quite low reciprocal dispersion (9A/mm) of their plate mate- 
rial, which is possibly not sufficient to perform high accuracy 
splitting measurements of the spectral lines observed in opposite 
polarization; ii) the employment of the Nasmyth focus station for 
the polarimeter: the oblique reflection from the tertiary mirror is 
prone to introduce a spurious polarization signal. This problem 
is much less critical at the Cassegrain focus, which was used in- 
stead for both the MuSiCoS observations obtained by Ch adid et] 
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[in](|2b04 ), and the FORS 1 measurements presented in this work. 
We therefore suggest that the magnetic field in RR Lyr was sig- 
nificantly overestimated in the earlier studies. 

5. Conclusion 

The results of our survey of magnetic fields in RR Lyrae stars 
reveal a serious challenge to the magnetic models for explain- 
ing the Blazhko effect. For a sample consisting of 17 RR Lyrae 
stars with different pulsational properties, no substantial dipole 
magnetic field, as required by the magnetic model (Shibahashi 
2000| l, was detected. We determined an upper limit for the 
strength of the dipole field in the stars of our sample at 130 G 
with a 95% confidence level. Our result implies that the Blazhko 
modulation in the pulsation of RR Lyrae stars is not correlated 
with the presence of a strong, quasi-dipolar magnetic field. For 
morphologies of higher order than the dipole field, however, 
our determinations of the mean longitudinal field as obtained 
in this work only set fairly loose constraints on the magnetic 
field strength. More complex magnetic field morphologies may 
be detected with high-resolution spectropolarimetric data. 
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